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ABSTRACT 

Analytical methods for the development of Reynolds stress models in turbulence are 
reviewed in detail. Zero, one and two equation models are discussed along with second-order 
closures. A strong case is made for the superior predictive capabilities of second-order closure 
models in comparison to the simpler models. The central points of the paper are illustrated 
by examples from both homogeneous and inhomogeneous turbulence. A discussion of the 
author’s views concerning the progress made in Reynolds stress modeling is also provided 
along with a brief history of the subject. 
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INTRODUCTION 


Despite over a century of research, turbulence remains the major unsolved problem of 
classical physics. While most researchers agree that the essential physics of turbulent flows 
can be described by the Navier-Stokes equations, limitations in computer capacity make 
it impossible - for now and the foreseeable future - to directly solve these equations in 
the complex turbulent flows of technological interest. Hence, virtually all scientific and 
engineering calculations of non-trivial turbulent flows, at high Reynolds numbers, are based 
on some type of modeling. This modeling can take a variety of forms: 


(a) Reynolds stress models which allow for the calculation of one-point first and second 
moments such as the mean velocity, mean pressure and turbulent kinetic energy. 

(b) Subgrid-scale models for large-eddy simulations wherein the large, energy containing 
eddies are computed directly and the effect of the small scales - which are more 
universal in character - are modeled. 

(c) Two-point closures or spectral models which provide more detailed information about 
the turbulence structure since they are based on the two-point velocity correlation 
tensor. 

(d) Pdf models based on the joint probability density function. 


Large-eddy simulations (LES) have found a variety of important geophysical applications 
where they have been used in weather forecasting as well as in other atmospheric studies (cf. 
Deardorff 1973, Clark and Farley 1984, and Smolarkiewicz and Clark 1985). Likewise, LES 
has shed new light on the physics of certain basic turbulent flows - which include homoge- 
neous shear flow and channel flow - at higher Reynolds numbers that are not accessible to di- 
rect simulations (cf. Moin and Kim 1982, Bardina, Ferziger and Reynolds 1983, Rogallo and 
Moin 1984, and Piomelli, Ferziger and Moin 1987). Two-point closures such as the EDQNM 
model of Orszag (1970) have been quite useful in the analysis of homogeneous turbulent 
flows where they have provided new information on the structure of isotropic turbulence (cf. 
Lesieur 1987) and on the effect of shear and rotation (cf. Bertoglio 1982). However, there 
are a variety of theoretical and operational problems with two-point closures and large-eddy 
simulations that make their application to strongly inhomogeneous turbulent flows difficult, 
if not impossible - especially in irregular geometries with solid boundaries. There have been 
no applications of two-point closures to wall-bounded turbulent flows and virtually all such 
applications of LES have been in simple geometries where Van Driest damping could be used 
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- an empirical approach that generally does not work well when there is flow separation. 
Comparable problems in dealing with wall-bounded flows have, for the most part, limited 
pdf methods to free turbulent flows where they have been quite useful in the description 
of chemically reacting turbulence (see Pope 1985). Since most practical engineering flows 
involve complex geometries with solid boundaries - at Reynolds numbers that are far higher 
than those that are accessible to direct simulations - the preferred approach has been to base 
such calculations on Reynolds stress modeling.* This forms the motivation for the present 
review paper whose purpose is to put into perspective some of the more recent theoretical 
developments in Reynolds stress modeling. 

The concept of Reynolds averaging was introduced by Sir Osborne Reynolds in his land- 
mark turbulence research of the latter part of the nineteenth century (see Reynolds 1895). 
During a comparable time frame, Boussinesq (1877) introduced the concept of the turbulent 
or eddy viscosity as the basis for a simple time-averaged turbulence closure. However, it 
was not until after 1920 that the first successful calculation of a practical turbulent flow was 
achieved based on the Reynolds averaged Navier-Stokes equations with an eddy viscosity 
model. This was largely due to the pioneering work of Prandtl (1925) who introduced the 
concept of the mixing length as a basis for the determination of the eddy viscosity. This mix- 
ing length model led to closed form solutions for turbulent pipe and channel flows that were 
remarkably successful in collapsing the existing experimental data. A variety of turbulence 
researchers — most notably including Von Karmdn (1930, 1948) — made further contributions 
to the mixing length approach which continued to be a highly active area of research until 
the post World War II period. By this time it was clear that the basic assumptions behind 
the mixing length approach - which makes a direct analogy between turbulent transport 
processes and molecular transport processes - were unrealistic; turbulent flows do not have 
a clear cut separation of scales. With the desire to develop more general models, Prandtl 
and Wieghardt (1945) tied the eddy viscosity to the turbulent kinetic energy which was 
obtained from a separate modeled transport equation. This was a precursor to the one equa- 
tion models of turbulence - or so called K — l models - wherein the turbulent length scale 
i is specified empirically and the turbulent kinetic energy K is obtained from a modeled 
transport equation. However, these models still suffered from the deficiencies intrinsic to all 
eddy viscosity models: the inability to properly account for streamline curvature and history 
effects on the individual Reynolds stress components. 

In a landmark paper by Rotta (1951), the foundation was laid for a full Reynolds stress 
turbulence closure which was to ultimately change the course of Reynolds stress modeling. 

t In fact, the only alternative of comparable simplicity is the vorticity transport theory of Taylor (1915); 
a three-dimensional vorticity covariance closure along these lines has been recently pursued by Bernard and 
co-workers (cf. Bernard and Berger 1982), 
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This new approach of Rotta - which is now referred to as second order or second moment 
closure - was based on the Reynolds stress transport equation. By making use of some of the 
statistical ideas of Kolmogorov from the 1940’s - and by introducing some entirely new ideas 
- Rotta succeeded in closing the Reynolds stress transport equation. This new Reynolds 
stress closure, unlike eddy viscosity models, accounted for both history and nonlocal effects 
on the evolution of the Reynolds stress tensor - features whose importance had long been 
known. However, since this approach required the solution of an additional six transport 
equations for the individual components of the Reynolds stress tensor, it was not to be 
computationally feasible for the next few decades to solve complex engineering flows based 
on a full second-order closure. By the 1970’s, with the wide availability of high speed 
computers, a new thrust in the development and implementation of second-order closure 
models began with the work of Daly and Harlow (1970) and Donaldson (1972). In an 
important paper, Launder, Reece and Rodi (1975) developed a new second-order closure 
model that improved significantly on the earlier work of Rotta (1951). This paper developed 
more systematic models for the pressure-strain correlation and turbulent transport terms; a 
modeled transport equation for the turbulent dissipation rate was solved in conjunction with 
this Reynolds stress closure. However, more importantly, Launder, Reece and Rodi (1975) 
showed how second-order closure models could be cahbrated and applied to the solution of 
practical engineering flows. When the Launder, Reece and Rodi (1975) model is contracted 
and supplemented with an eddy viscosity representation for the Reynolds stress, a two- 
equation model (referred to as the K — 6 model) is obtained which is identical to that 
derived by Hanjalic and Launder (1972) a few years earlier. Because of the substantially 
lower computational effort required, the K — e model is still one of the most commonly used 
turbulence models for the solution of practical engineering problems. 

Subsequent to the publication of the paper by Launder, Reece and Rodi (1975), a variety 
of turbulence modelers have continued research on second-order closures. Lumley (1978) in- 
troduced the important constraint of realizability and made significant contributions to the 
modeling of the pressure-strain correlation. Launder and co-workers continued to expand 
on the refinement and application of second-order closure models to problems of significant 
engineering interest (see Launder 1989). Speziale (1985, 1987a) exploited invariance argu- 
ments - along with consistency conditions for solutions of the Navier-Stokes equations in 
a rapidly rotating frame - to develop new models for the rapid pressure-strain correlation. 
Haworth and Pope (1986) developed a second-order closure model starting from the pdf 
based Langevin equation. Reynolds (1988) has attempted to develop models for the rapid 
pressure-strain correlation by using Rapid Distortion Theory. 

In this paper, analytical methods for the derivation of Reynolds stress models will be 
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reviewed. Zero, one and two-equation models will be considered along with second-order 
closures. Two approaches to the development of models will be discussed: 

(1) The continuum mechanics approach which is typically based on a Taylor expansion. 
Invariance constraints - as well other consistency conditions such as Rapid Distortion 
Theory (RDT) and realizability - are then used to simplify the model. The remaining 
constants are evaluated by reference to benchmark physical experiments. 

(2) The statistical mechanics approach which is based on the construction of an asymp- 
totic expansion. Unlike in the continuum mechanics approach, here the constants 
of the model are calculated explicitly. The two primary examples of this approach 
are the two-scale DIA models of Yoshizawa (1984) and the Renormalization Group 
(RNG) models of Yakhot and Orszag (1986). 


The basic methodology of these two techniques will be examined, however, more emphasis 
will be placed on the continuum mechanics approach since there is a larger body of literature 
on this method and since it has been the author’s preferred approach. The strengths and 
weaknesses of a variety of Reynolds stress models will be discussed in detail and illustrated 
by examples. A strong case will be made for the superior predictive capabilities of second- 
order closures in comparison to the older zero, one and two-equation models. However, 
some significant deficiencies in the structure of second-order closures that still remain will be 
pointed out. These issues, as well as the author’s views concerning possible future directions 
of research, will be discussed in the sections to follow. 


BASIC EQUATIONS OF REYNOLDS STRESS MODELING 


We will consider the turbulent flow of a viscous, incompressible fluid with constant prop- 
erties (limitations of space will not allow us to discuss compressible turbulence modeling in 
any detail). The governing field equations are the Navier-Stokes and continuity equations 
which are given by 

i • rli/ • Hn 

( 1 ) 


dui dm dp 2 

IT + u ’dTj = + " v Ui 


dui 

dxi 


= 0 


( 2 ) 


where is the velocity vector, p is the modified pressure (which can include a gravitational 
potential), and v is the kinematic viscosity of the fluid. In (1) - (2), the Einstein summation 
convention applies to repeated indices. 
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The velocity and pressure are decomposed into mean and fluctuating parts as follows: 

ui = u i + u' i , p = p + p'. (3) 

It is assumed that any flow variables <f> and ip obey the Reynolds averaging rules (cf. 
Tennekes and Lumley 1972): 

J' = ^ = 0 (4) 

= ~<p "ip + ffij? (5) 

<p'ip = ip'<p — 0. (6) 

In a statistically steady turbulence, the mean of a flow variable <f> can be taken to be the 
simple time average 

= ( 7 ) 

whereas for a spatially homogeneous flow, a volume average can be used 

1> = <? V) = y lim, p: J y ^(x, t)d 3 x. (8) 

For more general turbulent flows that are neither statistically steady nor homogeneous, the 
mean of any flow variable <p is taken to be the ensemble mean 

7 = ?'*>= lim W 

N ^°° N k= 1 

where an average is taken over N repeated experiments. The ergodic hypothesis is assumed 
to apply - namely, in a statistically steady turbulent flow it is assumed that 

-fi T) =$ E) (10) 


and in a homogeneous turbulent flow it is assumed that 

f V) = -$ E) . 


( 11 ) 


The Reynolds equation - which physically corresponds to a balance of mean linear mo- 
mentum - takes the form 


dui _ du{ 

~di + Uj d^ 


dp , y 7 2 — dT a 

~a~ + "V Ui- — 
OXi OXj 


( 12 ) 


where 



is the Reynolds stress tensor. Equation (12) is obtained by substituting the decompositions 
(3) into the Navier-Stokes equation (1) and then taking an ensemble mean. The mean 
continuity equation is given by 

0 (14) 


duj 

dxi 


and is obtained by simply taking the ensemble mean of (2). Equations (12) - (14) do not 
represent a closed system for the determination of the mean velocity Ui and mean pressure 
p due to the additional six unknowns contained within the Reynolds stress tensor. The 
problem of Reynolds stress closure is to tie the Reynolds stress tensor to the mean velocity 
field in some physically consistent fashion. 

In order to gain greater insight into the problem of Reynolds stress closure, we will 
now consider the governing field equations for the fluctuation dynamics. The fluctuating 
momentum equation - from which u\ is determined - takes the form 


M + iM _ u .p. _W + „ V “u; + dT,i 

Qj- 1 ft ? *_• 3 Bt.i 3 dx; dx; 


, .. T , (15) 

'dxj ~ '“iWj ~ " 3 dxj dxr ' dx d K J 

and is obtained by subtracting (12) from (1) after the decompositions (3) are introduced. 
The fluctuating continuity equation, which is obtained by subtracting (14) from (2), is given 

by 

(16) 


£- 0 . 

OXi 


Equations (15) - (16) have solutions for the fluctuating velocity u[ that are of the general 
mathematical form 


u'(x,t) = JPi[u(y,3),u , (y,0),u'(y,3)|av;x,t] y € V, s G (—«>,*) (17) 

where T x [ ■ } denotes a functional, V is the volume of the fluid, and dV is its bounding 
surface. In alternative terms, the fluctuating velocity is a functional of the global history 
of the mean velocity field with an implicit dependence on its own initial and boundary 
conditions. Here we use the term functional in its broadest mathematical sense, namely, any 
quantity determined by a function. From (17), we can explicitly calculate the Reynolds stress 
tensor ry = ujuj which will also be a functional of the global history of the mean velocity. 
However, there is a serious problem in regard to the dependence of on the initial and 
boundary conditions for the fluctuating velocity as discussed by Lumley (1970). There is no 
hope for a workable Reynolds stress closure if there is a detailed dependence on such initial 
and boundary conditions. For turbulent flows that are sufficiently far from solid boundaries 
- and sufficiently far evolved in time past their initiation - it is not unreasonable to assume 
that the initial and boundary conditions on the fluctuating velocity (beyond those for Tifi) 
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merely set the length and time scales of the turbulence. Hence, with this crucial assumption, 
we obtain the expression 


Tij(x^) = jF ii [u(y,s) ) 4(y,a),ro(y,5);x,<] y € V, s G (—oo,t) (18) 


where i 0 is the turbulent length scale, r 0 is the turbulent time scale, and the functional 
depends implicitly on the initial and boundary conditions for (see Lumley 1970 for a 
more detailed discussion of these points). Equation (18) serves as the cornerstone of Reynolds 
stress modeling. Eddy viscosity models, which are of the form 



(where the turbulent or eddy viscosity vt oc ^q/ t o) represent one of the simplest examples 
of (18). 

Since we will be discussing second-order closure models later, it would be useful at this 
point to introduce the Reynolds stress transport equation as well as the turbulent dissipation 
rate transport equation. The latter equation plays an important role in many commonly used 
Reynolds stress models where the turbulent dissipation rate is used to build up the turbulent 
length and time scales. If we denote the fluctuating momentum equation (15) in operator 
form as 

Cu\ = 0, (20) 

then the Reynolds stress transport equation is obtained from the second moment 


u'iCu'j + u'jCu'i = 0 


( 21 ) 


whereas the turbulent dissipation rate is obtained from the moment 


^ du' { d 
dxj dxj 


(£“!) = o. 


(22) 


More explicitly, the Reynolds stress transport equation (21) is given by (cf. Hinze 1975) 


dr. 


I TT 


dt 


+ Uk 


dTij duj &Ui dCijk , „ a 

dT t = ~ r *W k ~ T ’ k a7k + ” Cii “ aST + T ' 


*3 


(23) 


where 



_ du[ du'j 
V dxk dxk 


Cijk = u' l u' :j u' k + p'u’Jjk + 


(24) 

(25) 

(26) 
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are the pressure-strain correlation, dissipation rate correlation and third-order diffusion cor- 
relation, respectively. On the other hand, the turbulent dissipation rate transport equation 
(22) is given by 


de _ de 
dt U * dxi 


2 du[ du' k dui 

e ~ 2 v lh j dx, 


Jx^dv^d^ _ *4 d 2 u t 

dxi dx k dx k k dxj dx k dxj 

2 LMM) 

dxjc dx m dx m dxji y 9x m dx m J 


(27) 


_2„JL 

dx k dx m J dxkdx m dxkdxyyi 

where e = \eu is the scalar dissipation rate. The seven higher-order correlations on the 
right-hand-side of (27) correspond to three physical effects: the first four terms give rise 
to the production of dissipation, the next two terms represent the turbulent diffusion of 
dissipation, and the last term represents the turbulent dissipation of dissipation. 

Finally, before closing this section, it would be useful to briefly discuss two constraints 
that have played a central role in the formulation of modern Reynolds stress models: real- 
izability and frame invariance. The constraint of realizability was rigorously introduced by 
Lumley (see Lumley 1978, 1983 for a more detailed discussion). It requires that a Reynolds 
stress model yield positive component energies, i.e., that 

T*a>0, a = 1,2,3 (28) 

for any given turbulent flow. The inequality (28) (where Greek indices are used to indicate 
that there is no summation) is a direct consequence of the definition of the Reynolds stress 
tensor given by (13). It was first shown by Lumley that realizability could be satisfied 
identically in homogeneous turbulent flows by Reynolds stress transport models; this is 
accomplished by requiring that whenever a component energy T aa vanishes, its time rate T aa 
also vanishes. 

Donaldson (1968) was probably the first to advocate the unequivocal use of coordinate 
invariance in turbulence modeling. This approach, which Donaldson termed “invariant mod- 
eling,” was based on the Reynolds stress transport equation and required that all modeled 
terms be cast in tensor form. Prior to the 1970’s it was not uncommon for turbulence models 
to be proposed that were incapable of being uniquely put in tensor form (hence, these older 
models could not be properly extended to more complex flows, particularly to ones involv- 
ing curvilinear coordinates). The more complicated question of frame invariance — where 
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time-dependent rotations and translations of the reference frame are accounted for - was 
first considered by Lumley (1970) in an interesting paper. A more comprehensive analysis 
of the effect of a change of reference frame was conducted by the author in a series of papers 
published during the 1 980’s (see Speziale 1989a for a detailed review of these results). In an 
arbitrary non-inertial reference frame, which can undergo arbitrary time-dependent rotations 
and translations relative to an inertial framing, the fluctuating momentum equation takes 
the form 


du[ du[ . du' . dUi dp' 2 , dr^ , 

^ : + 4 " iik A 


at 


(29) 


where e tJ fc is the permutation tensor and is the rotation rate of the reference frame relative 
to an inertial framing (see Speziale 1989a). From (29), it is clear that the evolution of the 
fluctuating velocity only depends directly on the' motion of the reference frame through 
the Coriolis acceleration; translational accelerations - as well as centrifugal and angular 
accelerations - only have an indirect effect through the changes that they induce in the 
mean velocity field. Consequently, closure models for the Reynolds stress tensor must be 
form invariant under the extended Galilean group of transformations 


x* = x + c (t) 


(30) 


which allows for an arbitrary translational acceleration - c of the reference frame relative to 
an inertial framing x. 

In the limit of two-dimensional turbulence (or a turbulence where the ratio of the fluc- 
tuating to mean time scales tq/Tq < 1), the Coriolis acceleration is derivable from a scalar 
potential that can be absorbed into the fluctuating pressure (or neglected) yielding complete 
frame-indifference (see Speziale 1981, 1983). This invariance under arbitrary time-dependent 
rotations and translations of the reference frame specified by 


x* = Q(t)x + c(f) (31) 

(where Q(t) is any time-dependent proper- orthogonal rotation tensor) is referred to as Mate- 
rial Frame Indifference (MFI) - the term that has been traditionally used for the analogous 
manifest invariance of constitutive equations in continuum mechanics. For general three- 
dimensional turbulent flows where to/To = 0(1), MFI does not apply as first pointed out by 
Lumley (1970). However, the Coriolis acceleration in (29) can be combined with the mean 
velocity in such a way that frame-dependence enters exclusively through the appearance of 
the intrinsic or absolute mean vorticity defined by (see Speziale 1989a) 

Wij = 2 “ ft?) + (32) 
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This result, along with the constraint of MFI in the two-dimensional limit, restricts the 
allowable form of models considerably. 


ZERO-EQUATION AND ONE-EQUATION MODELS BASED ON AN EDDY 
VISCOSITY 


In the simplest continuum mechanics approach - whose earliest formulations have often 
been referred to as phenomenological models - the starting point is equation (18). Invariance 
under the extended Galilean group of transformations (30) - which any physically sound 
Reynolds stress model must obey - can be satisfied identically by models of the form 


Tijfat) = ^i[u(y,3)-Tr(x,s),4(y,s),r 0 (y,s);x,<] y e V, s6(-oo,f). (33) 


The variables u(y,s) — u(x, s 
follows: 

u( y,s) -H(x,s) 


),4>(y.a) an d T 0 (y, s) can be expanded in a Taylor series as 


% 9“ . (ft - ZiXft - x i) 9 3 u , 

- (vi “ Ii W + 2! a^ + 

3 2 u 


(34) 


.. N , , M , ^dlo (s-t) 2 dH 0 

lo{y,s) = lo + {y i -x i )Q- + (3 -t)—+ — — + 

+{s _ t){yi _ Xi) 9% + 


2 ! 


dxidxj 


dtdxi 


(35) 


7b(y, a) = r 0 + (yi - Xi) + (a - i)“ + 

(Vi - Xi)(yj - Xj) d 2 r 0 


dr 0 ( 3 — t) 2 d 2 r 0 


2! di 2 


+ 


pjl 

+ (3- t)( Vi - + 


(36) 


2! dx^ ' 'dtdxi 

where terms up to the second order are shown and it is understood that U, £ 0 and r 0 on the 
r.h.s. of (34) - (36) are evaluated at x and t. After splitting into isotropic and deviatoric 
parts - and applying elementary dimensional analysis - the following expression is obtained: 


2 l 2 

Tij = -KSij - -|£j[v(y, 3 ) - v(x, s)\ x, t] y G V, (- 00 , t) 


(37) 


where 


_ r 0 u 

v = ~r> 

to 


K = -TV, 


(38) 
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are, respectively, the dimensionless mean velocity and the turbulent kinetic energy. Tij is 
a traceless and dimensionless functional of its arguments. By making use of the Taylor 
expansions (34) - (36), it is a simple matter to show that 


-to 


where 


* * _ W ~ X < 

Vi ~ x i = 


to 


! suj y agi 

\dxj ) °dzj 


(39) 


(40) 


are dimensionless variables of order one given that To is the time scale of the mean flow. If, 
analogous to the molecular fluctuations of most continuum flows, we assume that there is a 
complete separation of scales such that 

to 


-to 


Lq 


< 1, 


(41) 


equation (37) can then be localized in space and time. Of course, it is well-known that this 
constitutes an over-simplification; the molecular fluctuations of most continuum flows are 
such that Tq/Tq < 10 -6 whereas with turbulent fluctuations, r 0 /T 0 can be of 0(1). 

By making use of (39) - (41), equation (37) can be localized to the approximate form 


r ii = \K6ij - 4g«( v kf ) 

O Tq 


where 


Vk,t = 


T 0 


'duu 


(42) 


(43) 


T 0 \dx t J 

is the dimensionless mean velocity gradient. Since the tensor function Gij is symmetric 
and traceless (and since Vij is traceless) it follows that — to the first order in tq/Tq - form 
invariance under a change of coordinates simplifies (42) to (cf. Smith 1971): 


2 (dui duA 


(44) 


where 

v t=^o/ t o (45) 

is the eddy viscosity. While the standard eddy viscosity model (44) comes out of this 
derivation when only first order terms in To/To are maintained, anisotropic eddy viscosity 
(or viscoelastic) models are obtained when second-order terms axe maintained. These more 
complicated models will be discussed in the next section. 

Eddy viscosity models are not closed until prescriptions are made for the turbulent length 
and time scales in (45). In zero equation models, both Iq and r 0 are prescribed algebraically. 
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The earliest example of a successful zero-equation model is PrandtPs mixing length theory 
(see Prandtl 1925). By making analogies between the turbulent length scale and the mean 
free path in the kinetic theory of gases, Prandtl argued that vt should be of the form 


Vj 


e |*| 

ml dy ] 


(46) 


for a plane shear flow where the mean velocity is of the form u = tt(y) i. In (46), £ m is 
the “mixing length” which represents the distance traversed by a small lump of fluid before 
losing its momentum. Near a plane solid boundary, it was furthermore assumed that 


lm = *y 


(47) 


where k is the Von Karm&n constant (this result can be obtained from a first-order Taylor 
series expansion since l m must vanish at a wall). When (46) - (47) are used in conjunction 
with the added assumption that the shear stress is approximately constant in the near wall 
region, the celebrated “law of the wall” is obtained: 


u + = —in y + + C 


K 


(48) 


where y + is measured normal from the wall and 


+ u , yu r 

« = — , y = — 

t I t V 


(49) 


given that u T is the friction velocity and C is a dimensionless constant. Equation (48) (with 
k = 0.4 and C = 4.9) was remarkably successful in collapsing the experimental data for 
turbulent pipe and channel flows for a significant range of t/ + varying from 30 to 1,000 (see 
Schlichting 1968 for an interesting review of these results). The law of the wall is still heavily 
used to this day as a boundary condition in the more sophisticated turbulence models for 
which it is either difficult or too expensive to integrate directly to a solid boundary. 

During the 1960’s and 1970’s, with the dramatic emergence of computational fluid dy- 
namics, some efforts were made to generalize mixing length models to three-dimensional tur- 
bulent flows. With such models, Reynolds averaged computations could be conducted with 
any existing Navier- Stokes computer code that allowed for a variable viscosity. PrandtPs 
mixing length theory (46) has two straightforward extensions to three dimensional flows: 
the strain rate form 

VT = o (50) 


where ~Sij = \(du i/dxj + duj/dxi) is the mean rate of strain tensor, or the vorticity form 


vt 



(51) 
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where &i = e^du^/ dxj is the mean vorticity vector. The former model (50) is due to 
Smagorinsky (1963) and has been primarily used as a subgrid scale model for large-eddy 
simulations; the latter model (51) is due to Baldwin and Lomax (1978) and has been widely 
used for Reynolds averaged aerodynamic computations. Both models - which collapse to 
Prandtl’s mixing length theory (46) in a plane shear flow - have the primary advantage of 
their computational ease of application. They suffer from the disadvantage of the need for an 
ad hoc prescription of the turbulent length scale in each problem solved as well as from the 
complete neglect of history effects. Furthermore, they do not provide for the computation 
of the turbulent kinetic energy which is a crucial measure of the intensity of the turbulence 
(such zero-equation models only allow for the calculation of the mean velocity and mean 
pressure). 

One-equation models were developed in order to eliminate some of the deficiencies cited 
above, namely, to provide for the computation of the turbulent kinetic energy and to account 
for some limited nonlocal and history effects in the determination of the eddy viscosity. In 
these one-equation models of turbulence, the eddy viscosity is assumed to be of the form 
(see Kolmogorov 1942 and Prandtl and Wieghardt 1945) 

v T = Kh (52) 


where the turbulent kinetic energy K is obtained from a modeled version of its exact trans- 
port equation 


dK _ dK duj __d_ 

1 W + Ui dxi ~ _T,i ~ e ~ dxi 


(\ u k u k u i+P' u 'i) 


+ uV 2 K. 


(53) 


Equation (53), which is obtained by a simple contraction of (23), can be closed once models 
for the turbulent transport and dissipation terms (i.e., the second and third terms on the r.h.s. 
of (53)) are provided. Consistent with the assumption that there is a clear-cut separation of 
scales (i.e., that the turbulent transport processes parallel the molecular ones), the turbulent 
transport term is modeled by a gradient transport hypothesis, i.e., 


-u' k u' k u( + p'u' x 


vt dK 

a K dxi 


(54) 


where <jk is a dimensionless constant. By simple scaling arguments - analogous to those 
made by Kolmogorov (1942) - the turbulent dissipation rate e is usually modeled as follows 


£ = C* 


K t 


(55) 


where C* is a dimensionless constant. A closed system of equations for the determination of 
Ui,p and K is obtained once the turbulent length scale i is specified empirically. It should 
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be mentioned that the modeled transport equation for the turbulent kinetic energy specified 
by equations (53) - (55) cannot be integrated to a solid boundary. Either wall functions 
must be used or low- Reynolds- number versions of (53) - (55) must be substituted (cf. Norris 
and Reynolds 1975 and Reynolds 1976). It is interesting to note that Bradshaw, Ferriss and 
Atwell (1967) considered an alternative one-equation model, based on a modeled transport 
equation for the Reynolds shear stress uV, which seemed to be better suited for turbulent 
boundary layers. 

Since zero and one equation models have not been in the forefront of turbulence modeling 
research for the past twenty years, we will not present the results of any illustrative calcu- 
lations (the reader is referred to Cebeci and Smith 1974 and Rodi 1980 for some interesting 
examples). The primary deficiencies of these models are twofold: (a) the use of an eddy 
viscosity, and (b) the need to provide an ad hoc specification of the turbulence length scale. 
This latter deficiency in regard to the length scale makes zero and one equation models 
incomplete; the two-equation models that will be discussed in the next section were the first 
complete turbulence models (i.e., models that only require the specification of initial and 
boundary conditions for the solution of problems). Nonetheless, despite these deficiencies, 
zero and one equation models have made some important contributions to the computation 
of practical engineering flows. Their simplicity of structure — and reduced computing times - 
continue to make them the most commonly adopted models for complex aerodynamic calcu- 
lations (see Cebeci and Smith 1968 and Johnson and King 1984 for two of the most popular 
such models). 

TWO-EQUATION MODELS 

A variety of two- equation models — which are among the most popular Reynolds stress 
models for scientific and engineering calculations — will be discussed in this section. Models of 
the K — e, K — l and K — u type will be considered based on an isotropic and anisotropic eddy 
viscosity. Both the continuum mechanics and statistical mechanics approach for deriving 
such two-equation models will be discussed. 

The feature that distinguishes two-equation models from zero- or one-equation models is 
that two separate modeled transport equations are solved for the turbulent length and time 
scales (or for any two linearly independent combinations thereof). In the standard K — e 
model — which is probably the most popular such model — the length and time scales are 
built up from the turbulent kinetic energy and dissipation rate as follows (see Hanjalic and 
Launder 1972 and Launder and Spalding 1974): 

, K\ K 

£q OC , To oc — . 

£ £ 
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Separate modeled transport equations are solved for the turbulent kinetic energy K and 
turbulent dissipation rate e. In order to close the exact transport equation for K, only a 
model for the turbulent transport term on the r.h.s. of (53) is needed; consistent with the 
overriding assumption that there is a clear-cut separation of scales, the gradient transport 
model (54) is used. The exact transport equation for e, given by (27), can be rewritten in 
the form 

' = uV 2 e + V e -1- V e - (56) 


de de 


where V e represents the production of dissipation (given by the first four correlations on the 
r.h.s. of (27)), V e represents the turbulent diffusion of dissipation (given by the next two 
correlations on the r.h.s. of (27)), and represents the turbulent dissipation of dissipa- 
tion (given by the last term on the r.h.s. of (27)). Again, consistent with the underlying 
assumption (41), a gradient transport hypothesis is used to model V e : 


v. = — ( 

dxi \ CT e dxi ) 


(57) 


where a e is a dimensionless constant. The production of dissipation and dissipation of 
dissipation are modeled as follows: 


pfu- 

V e = V e (b tj ,^-,K,s) 


(58) 

*. = $.(JT,c). (59) 

Eqs. (58) - (59) are based on the physical reasoning that the production of dissipation is 
governed by the level of anisotropy b±j = (r tJ - — &K6ij)/2K in the Reynolds stress tensor and 
the mean velocity gradients (scaled by K and e which determine the length and time scales) 
whereas the dissipation of dissipation is determined by the length and time scales alone (an 
assumption motivated by isotropic turbulence). By a simple dimensional analysis it follows 
that 2 

$ e = Cj- (60) 

where (7*2 is a dimensionless constant. Coordinate invariance coupled with a simple dimen- 
sional analysis yields 


V e - 2C el ebij 


(61) 


= -C el ^T, 


e du{ 

K Tij d7j 


as the leading term in a Taylor expansion of (58) assuming that ||b|| and tq/T q are small 
(C e i is a dimensionless constant). Equation (61) was originally postulated based on the 
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simple physical reasoning that the production of dissipation should be proportional to the 
production of turbulent kinetic energy (cf. Hanjalic and Launder 1972). A composition of 
these various modeled terms yields the standard K — £ model (cf. Launder and Spalding 
1974): 



Here, the constants assume the approximate values of C ^ = 0.09, &k = 1-0, cr e = 1.3, C e \ = 
1.44 and C e2 = 1-92 which are obtained, for the most part, by comparisons of the model 
predictions with the results of physical experiments on equilibrium turbulent boundary layers 
and the decay of isotropic turbulence. It should be noted that the standard K — £ model 
(62) cannot be integrated to a solid boundary; either wall functions or some form of damping 
must be implemented (see Patel, Rodi and Scheuerer 1985 for an extensive review of these 
methods). 

At this point, it would be useful to provide some examples of the performance of the 
K — s model in some benchmark, homogeneous turbulent flows as well as in a non-trivial, 
inhomogeneous turbulent flow. It is a simple matter to show that in isotropic turbulence 
where 

T ij = ~ £ij = — e(t)8ij 

the K — e model predicts the following rate of decay of the turbulent kinetic energy (cf. 
Reynolds 1987): 

K(t) = JST 0 [1 + (<7,2 - l)e 0 (63) 

Equation (63) indicates a power law decay where K ~ t~ 1-1 - a result that is not far removed 
from what is observed in physical experiments (cf. Comte-Bellot and Corrsin 1971). 

Homogeneous shear flow constitutes another classical turbulent flow that has been widely 
used to evaluate models. In this flow, an initially isotropic turbulence is subjected to a 
constant shear rate S with mean velocity gradients 


dui 

dxj 


0 5 0 \ 
0 0 0 
0 0 0 / 
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(64) 



The time evolution of the turbulent kinetic energy obtained from the standard K — e model 
is compared in Figure 1 with the large-eddy simulation of Bardina, Ferziger and Reynolds 
(1983) (here, K* = Kj Kq\s the dimensionless kinetic energy and t* = St is the dimensionless 
time). In so far as the equilibrium states are concerned, the standard K — e model predicts 
that (see Speziale and Mac Giolla Mhuiris 1989a): (6x2)00 = —0.217 and ( SK/e )«, = 4.82 in 
comparison to the experimental values of (612)00 = —0.15 and (SK/e)oo = 6.08, respectively* 
(see Tavoularis and Corrsin 1981). Consistent with a wide range of physical and numerical 
experiments, the standard K — e model predicts that the equilibrium structure of homoge- 
neous shear flow is universal (i.e., attracts all initial conditions) in the phase space of 6,y and 
SK/e. Hence, from Figure 1 and the equilibrium results given above, it is clear that the 
K — e model yields a qualitatively good description of shear flow; the specific quantitative 
predictions, however, could be improved upon. 

As an example of the performance of the standard K — e model in a more complicated 
inhomogeneous turbulence, the case of turbulent flow past a backward facing step at a 
Reynolds number Re ~ 100,000 will be presented (the same test case considered at the 
1980/81 AFOSR-HTTM Stanford Conference on Turbulence; it corresponds to the experi- 
mental test case of Kim, Kline, and Johnston 1980). In Figures 2(a) - (b) the mean flow 
streamlines and turbulence intensity profiles predicted by the K — e model are compared 
with the experimental data of Kim, Kline, and Johnston (1980). The standard K — e model 
- integrated using a single log wall layer starting at y + — 30 - predicts a reattachment point 
of x/AH == 5.7 in comparison to the experimental mean value of xj AH = 7.0. This error, 
which is of the order of 20%, is comparable to that which occurs in the predicted turbu- 
lence intensities (see Figure 2(b) and Speziale and Ngo 1988 for more detailed comparisons). 
However, Avva, Kline and Ferziger (1988) reported an improved prediction of x/AH = 6.3 
for the reattachment point by using a fine mesh and a three-layer log wall region. 

Recently, Yakhot and Orszag (1986) derived a K — e model based on RNG methods. 
In this approach, an expansion is made about an equilibrium state with known Gaussian 
statistics by making use of the correspondence principle. Bands of high wavenumbers (i.e., 
small scales) are systematically removed and space is rescaled. The dynamical equations for 
the renormalized (large-scale) velocity field account for the effect of the small scales that have 
been removed through the presence of an eddy viscosity. The removal of only the smallest 
scales gives rise to subgrid scale models for large-eddy simulations; the removal of successively 
larger and larger scales ultimately gives rise to Reynolds stress models. In the high Reynolds 
number limit, the RNG based K — e model of Yakhot and Orszag (1986) is identical in 
form to the standard K — e model (62). However, the constants of the model are calculated 

*Here, (-)oo denotes the equilibrium value obtained in the limit as t — » 00. 
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explicitly by the theory to be: = 0.0837,(7^ = 1.063, C e2 = 1.7215, <7^ = 0.7179 and 

cr, = 0.7179. Beyond having the attractive feature of no undetermined constants, the RNG 
K — e model of Yakhot and Orszag (1986) automatically bridges the eddy viscosity to the 
molecular viscosity as a solid boundary is approached, eliminating the need for the use of 
empirical wall functions or Van Driest damping. It must be mentioned, however, that some 
problems with the specific numerical values of the constants in the RNG K — e model have 
recently surfaced. In particular, the value of C c \ = 1.063 is dangerously close to C e \ = 1 
which constitutes a singular point of the e-transport equation. For example, the growth rate 
A of the turbulent kinetic energy (where K ~ e At * for \t* >■ 1) predicted by the K — e model 
in homogeneous shear flow is given by (see Speziale and Mac Giolla Mhuiris 1989a): 

, f c^C'2-c^y l* 

1(C,1 - 1 )(C« - 1). 

which becomes singular when C e \ = 1. Consequently, the value of C e i = 1.063 derived 
by Yakhot and Orszag (1986) yields excessively large growth rates for the turbulent kinetic 
energy in homogeneous shear flow in comparison to both physical and numerical experiments 
(see Speziale, Gatski and Mac Giolla Mhuiris 1989). 

One of the major deficiencies of the standard K —e model lies in its use of an eddy viscosity 
model for the Reynolds stress tensor. Eddy viscosity models have two major deficiencies 
associated with them: (a) they are purely dissipative and, hence, cannot account for Reynolds 
stress relaxation effects, and (b) they are oblivious to the presence of rotational strains (e.g., 
they fail to distinguish between the physically distinct cases of plane shear, plane strain, and 
rotating plane shear). In an effort to overcome these deficiencies, a considerable research 
effort has been directed toward the development of nonlinear or anisotropic generalizations of 
eddy viscosity models. By keeping second-order terms in the Taylor expansions (34) - (36), 
subject to invariance under the extended Galilean group (30), a more general representation 
for the Reynolds stress tensor is obtained: 



T a 


where 


\K6ii - 2 ^ + aJl (5 tk S kj - 

+ a 2 ^o jk& t k) 

1 / du{ dilj\ _ 1 / dui duA 

2 V&Cj dx{ ) ’ 2 \dxj dx< ) 


( 66 ) 


(67) 
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are the mean rate of strain and mean vorticity tensors (oq . . . cl 4 are dimensionless constants; 
in the linear limit as ol{ — > 0, the eddy viscosity model (44) is recovered). When a 4 = 0, the 
deviatoric part of (66) is of the general form o T ij = AijkiduiJ dxi (where Aijki depends al- 
gebraically on the mean velocity gradients) and, hence, the term “anisotropic eddy viscosity 
model” has been used in the literature. These models are probably more accurately char- 
acterized as “nonlinear” or “viscoelastic” corrections to the eddy viscosity models. Lumley 
(1970) was probably the first to systematically develop such models (with a 4 = 0) wherein he 
built up the length and time scales from the turbulent kinetic energy, turbulent dissipation 
rate, and the invariants of 15 ij and Z Jij. Saffman (1977) proposed similar anisotropic models 
which were solved in conjunction with modeled transport equations for K and to 2 (where 
u) = c/K). Pope (1975) and Rodi (1976) developed alternative anisotropic eddy viscosity 
models from the Reynolds stress transport equation by making an equilibrium hypothesis. 
Yoshizawa (1984, 1987) derived a more complete two-equation model - with a nonlinear 
correction to the eddy viscosity of the full form of (66) — by means of a two- scale Direct In- 
teraction Approximation (DIA) method. In this approach, Kraichnan’s DIA formalism (cf. 
Kraichnan 1964) is combined with a scale expansion technique where the slow variations of 
the mean fields are distinguished from the fast variations of the fluctuating fields by means 
of a scale parameter. The length and time scales of the turbulence are built up from the 
turbulent kinetic energy and dissipation rate for which modeled transport equations are de- 
rived. These modeled transport equations are identical in form to (62c) and (62d) except for 
the addition of higher-order cross diffusion terms. The numerical values of the constants are 
derived directly from the theory (as with the RNG K — s model). However, in applications 
it has been found that these values need to be adjusted (see Nisizima and Yoshizawa 1987). 

Speziale (1987b) developed a nonlinear K — e model based on a simplified version of (66) 
obtained by invoking the constraint of MFI in the limit of two-dimensional turbulence. In 
this model — where the length and time scales are built up from the turbulent kinetic energy 
and dissipation rate - the Reynolds stress tensor is modeled as^ 


'*7 



( 68 ) 


where 



dSij 

dt 


+ Uk 


dSij 

dxk 


dui &lij -g 

ox k dxk 


(69) 


§It is interesting to note that Rubinstein and Barton (1989) recently derived an alternative version of this 
model - which neglects the convective derivative in (69) -by using the RNG method of Yakhot and Orszag. 
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is the frame-indifferent Oldroyd derivative of "Sy and Cd — Ce = 1.68. Equation (68) can 
also be thought of as an approximation for turbulent flows where To /To -C 1 since MFI 
(which (68) satisfies identically) becomes exact in the limit as t 0 /T 0 — > 0. This model was 
shown by Speziale (1987b) and Speziale and Mac Giolla Mhuiris (1989a) to yield much 
more accurate predictions for the normal Reynolds stress anisotropies in turbulent channel 
flow and homogeneous shear flow (the standard K — e model erroneously predicts that 
r xx = Tyy = t zz = \K). As a result of this feature, the nonlinear K — e model is capable 
of predicting turbulent secondary flows in non- circular ducts unlike the standard K — e 
model which erroneously predicts a unidirectional mean turbulent flow (see Figures 3(a) - 
(c)). Comparably good predictions of turbulent secondary flows in a rectangular duct were 
obtained much earlier by Launder and Ying (1972), Gessner and Po (1976) and Demuren and 
Rodi (1984) who used the nonlinear algebraic Reynolds stress model of Rodi (1976). Due to 
the more accurate prediction of normal Reynolds stress anisotropies - and the incorporation 
of weak relaxation effects - the nonlinear K — e model of Speziale (1987b) was also able to 
yield improved results for turbulent flow past a backward facing step (compare Figures 4(a) - 
(b) with Figures 2(a) - (b)). Most notably, the nonlinear K — e model predicts reattachment 
at xj AH = 6.4 - a value which is more in line with the experimental value of xj AH = 7.0 
(as shown earlier, the standard K — e model yields a value of xj AH = 5.7 when a single log 
wall layer is used). 

Alternative two-equation models based on the solution of a modeled transport equation 
for an integral length scale (the K — l model) or the reciprocal time scale (the K — w model) 
have also been considered during the past fifteen years. In the K — l model (see Mellor and 
Herring 1973) a modeled transport equation is solved for the integral length scale i defined 

by 



-Rt«(x,r, 
47r |r| 2 



(70) 


where Rij = u((x, t)u'(x + r, t) is the two-point velocity correlation tensor. The typical form 
of the modeled transport equation for l is as follows: 


d(Kl) _ d(Kl) 

-(- Tij" 


dt 


dxi 


_d_ 

dxi 




(71) 


where /? x ,.../? 4 are empirical constants. Equation (71) is derived by integrating the con- 
tracted form of a modeled transport equation for the two-point velocity correlation tensor 
Ra (see Wolfshtein 1970). Mellor and co-workers have utilized this K — l model - with an 
eddy viscosity of the form (52) - in the solution of a variety of engineering and geophysical 
fluid dynamics problems (see Mellor and Herring 1973 and Mellor and Yamada 1974 for a 
more thorough review). It has been argued - and correctly so - that it is more sound to 
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base the turbulent macroscale on the integral length scale (70) rather than on the turbulent 
dissipation rate which only formally determines the turbulent microscale. However, for ho- 
mogeneous turbulent flows, it is a simple matter to show that this K —l model is equivalent 
to a K — e model where the constants C^, C e \ and C e 2 assume slightly different values (cf. 
Speziale 1989b). Furthermore, the modeled transport equation (71) for l has comparable 
problems to the modeled e-transport equation in so fax as integrations to a solid boundary 
are concerned (either wall functions or wall damping must be used). Consequently, at their 
current stage of development, it does not appear that this type of K — l model offers any 
significant advantages over the K — e model. 

Wilcox and co-workers have developed two-equation models of the K — lu type (see Wilcox 
and Traci 1976 and Wilcox 1988). In this approach, modeled transport equations are solved 
for the turbulent kinetic energy K and reciprocal turbulent time scale u> = e/K. The 
modeled transport equation for u> is of the form 


du> 

~dt 


, _ d(v u dui 2 d ( vt du> \ 

+ Ui d7i = ~ 7l K Tii d^ + !hi J + " v 


UJ 


(72) 


where v? = ~f*K/u and 71,72,7 * and cr u are constants. Equation (72) is obtained by 
making the same type of assumptions in the modeling of the exact transport equation for 
u that were made in developing the modeled e-transport equation (62d). For homogeneous 
turbulent flows, there is little difference between the K — e and K — oj models. The primary 
difference between the two models is in their treatment of the transport terms: the K — e 
model is based on a gradient transport hypothesis for e whereas the K — uj model uses the 
same hypothesis for oj instead. Despite the fact that ui is singular at a solid boundary, there 
is some evidence to suggest that the K — u> model is more computationally robust for the 
integration of turbulence models to a wall (i.e., there is the need for less empirical damping; 
see Wilcox 1988). 


SECOND-ORDER CLOSURE MODELS 

Theoretical Background 

Although two-equation models represent the first simple and complete Reynolds stress 
models to be developed, they still have significant deficiencies that make their application 
to complex turbulent flows precarious. As mentioned earlier, the two-equation models of 
the eddy viscosity type have the following major deficiencies: (a) the inability to properly 
account for streamline curvature, rotational strains and other body force effects and (b) the 
neglect of nonlocal and history effects on the Reynolds stress anisotropies. Most of these 
deficiencies are intimately tied to the assumption that there is a clear-cut separation of scales 


21 


at the second moment level (i.e., the level of the Reynolds stress tensor). This can be best 
illustrated by the example of homogeneous shear flow presented in the previous section. For 
this problem, the ratio of fluctuating to mean time scales is given by 



for the K — e model. This is in flagrant conflict with the assumption that t 0 /T 0 <C 1 which 
was crucial for the derivation of the K — £ model! While some of the deficiencies cited above 
can be partially overcome by the use of two-equation models with a nonlinear algebraic 
correction to the eddy viscosity, major improvements can only be achieved by higher-order 
closures - the simplest of which are second-order closure models. 

Second-order closure models are based on the Reynolds stress transport equation (23). 
Since this transport equation automatically accounts for the convection and diffusion of 
Reynolds stresses, second-order closure models (unlike eddy viscosity models) are able to 
account for strong nonlocal and history effects. Furthermore, since the Reynolds stress 
transport equation contains convection and production terms that adjust themselves auto- 
matically in turbulent flows with streamline curvature or a system rotation (through the 
addition of scale factors or Coriolis terms), complex turbulent flows involving these effects 
are usually better described. 

In order to develop a second-order closure, models must be provided for the higher-order 
correlations C ijk) II#, and e i3 on the right-hand-side of the Reynolds stress transport equation 
(23). These terms, sufficiently far from solid boundaries, are typically modeled as follows: 

(1) The third-order transport term C ijk is modeled by a gradient transport hypothesis 
which is based on the usual assumption that there is a clear-cut separation of scales 
between mean and fluctuating fields. 

(2) The pressure- strain correlation Iljy and the dissipation rate correlation are modeled 
based on ideas from homogeneous turbulence wherein the departures from isotropy 
are assumed to be small enough to allow for a Taylor series expansion about a state 
of isotropic turbulence. 

Near solid boundaries, either wall functions or wall damping are used in a comparable 
manner to that discussed m the last section. One important point to note is that the crucial 
assumption of separation of scales is made only at the third moment level. This leads us 
to the raison d’etre of second-order closure modeling: Since crude approximations for the 
second moments in eddy viscosity models can often yield adequate approximations for first 
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order moments (i.e., u and p) it may follow that crude approximations for the third-order 
moments can yield adequate approximations for the second-order moments in Reynolds stress 
transport models. 

The pressure-strain correlation 11^ plays a crucial role in determining the structure of 
most turbulent flows. Virtually all of the models for II jy that have been used in conjunction 
with second-order closure models are based on the assumption of local homogeneity. For 
homogeneous turbulent flows, the pressure-strain correlation takes the form 


where 


n»j A%j "h hdijki 


duk 

dxi 


(73) 


, i _ r <tv(y) 

1 4ir 7-oo \dxj dxi ) dy k dyi |x - y| 

_1_ r (dvj du‘Adu’ t dVM 
13 2n J-oo V&cy dxi ) dy k |x - y| ' 

Here, the first term on the r.h.s. of (73) is referred to as the slow pressure-strain whereas 
the second term is called the rapid pressure-strain. It has been shown that A^ and M tjki are 
functionals - in time and wavenumber space - of the energy spectrum tensor (cf. Weinstock 
1981 and Reynolds 1987). In a one-point closure, this suggests models for and Mjj kl 
that are functionals of the Reynolds stress tensor and turbulent dissipation rate. Neglecting 
history effects, the simplest such models are of the form 


Aij = eAjj( b) (74) 

Mijki = KMijkiih). (75) 

These algebraic models - based on the assumptions stated above - are obtained by using sim- 
ple dimensional arguments combined with the fact that II;y vanishes in the limit of isotropic 
turbulence (a constraint identically satisfied if ,4^(0) = 0 and Miju{ 0) = 0). Virtually all 
of the models for the pressure-strain correlation that have been used in conjunction with 
second-order closure models are of the form (73) - (75). 

Lumley (1978) was probably the first to systematically develop a general representation 
for the pressure-strain correlation based on (73) - (75). It can be shown that invariance under 
a change of coordinates - coupled with the assumption of analyticity about the isotropic state 
b{j = 0 - restricts (73) to be of the form^ 

^This representation, obtained by using the results of Smith (1971) on isotropic tensor functions, is 
actually somewhat more compact than that obtained by Lumley and co-workers. 
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Ily = aosbij + 0'\^{bikbkj — — 1 1 5 t j) + 0-2 K~Bij 

+(a 3 tr b-S + a 4 tr b 2 -S )Kbij + ( a s tr b-5 
+a 6 tr b 2 S )K(bikbkj — + °7 K{b ik Sj k 

+bjk^ik - \tT b-S Sij) + a s K(bikbkiSji 

O 

+bjkbki^u ~ g^ r b 2 >5 S{j) + ag K(bi k TDjk 


+bj k u ik ) + a w K(b lk b k iu}ji + b jk b k tu}u ) 


(76) 


where 

a; = Oj(//, ///), 1 = 0,1,... ,10 

II = bijbij, III — bikbkibii. 

and ir(-) denotes the trace. The eigenvalues b^ of bij are bounded as follows (see Lumley 
1978) 

- | < 6 ( “> < |, a = 1,2,3 (77) 

and for many engineering flow, ||b || 2 = |&^|max < 0.25. Hence, it would seem that a low 
order Taylor series truncation of (76) could possibly provide an adequate approximation. To 
the first order in b,ji 

n ij = -Ciebij-hC 2 KZij+C 3 K (bikSjk + bj&k - g bmn^ mn ) +C 4 K(b ik uJ jk +b jk U ik ) (78) 

which is the form used in the Launder, Reece and Rodi (1975) model. In the Launder, Reece 
and Rodi (LRR) model, the constants C\, C 3 , and C 4 were calibrated based on the results of 
return to isotropy and homogeneous shear flow experiments. The constant C 2 was chosen to 
be consistent with the value obtained by Crow (1968) from RDT for an irrotationally strained 
turbulence starting from an initially isotropic state. This yielded the following values for 
the constants in the simplified version of the LRR model: C\ = 3.6, C 2 = 0.8, C 3 — 0.6 and 
C4 = 0.6. It should be noted that the representation for the slow pressure-strain correlation 
in the LRR model is the Rotta (1951) return to isotropy model with the Rotta constant C\ 
adjusted from 2.8 to 3.6 (a value which is in the range of what can be extrapolated from 
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physical experiments). This model — consistent with experiments — predicts that an initially 
anisotropic, homogeneous turbulence relaxes gradually to an isotropic state after the mean 
velocity gradients are removed. 

An even simpler version of (78) was proposed by Rotta (1972) wherein <7 3 = (7 4 = 0. This 
model has been used by Mellor and co-workers for the calculation of a variety of engineering 
and geophysical flows (see Mellor and Herring 1973 and Mellor and Yamada 1974). Research 
during the past decade has focused attention on the development of nonlinear models for 
n iy Lumley (1978) and Shih and Lumley (1985) developed a nonlinear model by using the 
constraint of realizability discussed earlier. Haworth and Pope (1986) developed a nonlinear 
model for the pressure-strain correlation based on the Langevin equation used in the pdf 
description of turbulence. This model — which was cubic in the anisotropy tensor — was cali- 
brated based on homogeneous turbulence experiments and was shown to perform quite well 
for a range of such flows. Speziale (1987a) developed a hierarchy of second-order closure mod- 
els that were consistent with the MFI constraint in the limit of two-dimensional turbulence 
(this constraint was also made use of by Haworth and Pope 1986 in the development of 
their second- order closure). Launder and his co-workers (cf. Fu, Launder and Tselepidakis 
1987 and Craft, et al. 1989) have developed new nonlinear models for the pressure-strain 
correlation based on the use of realizability combined with a calibration using newer ho- 
mogeneous turbulence experiments. Reynolds (1988) has attempted to develop models that 
are consistent with RDT and the author has been analyzing models based on a dynamical 
systems approach (see Speziale and Mac Giolla Mhuiris 1989a, 1989b and Speziale, Sarkar 
and Gatski 1990). 

The modeling of the dissipation rate tensor, at high turbulence Reynolds numbers, is 
usually based on the Kolmogorov hypothesis of isotropy where 

Cij = (79) 


given that c == udu f ^j dxjdu^j dxj is the scalar dissipation rate. Here, the turbulent dissipation 
rate 6 is typically taken to be a solution of the modeled transport equation 


de _ dc 
dt U * dxi 


= -C el 


e dui 

K^’ThTj 


Cl2 lr +C, ^;( 


T a 



(80) 


where C el = 1.44, C e2 = 1.92 and C e = 0.15. Equation (80) is identical to the e-transport 
equation used in the K — e model with one exception: the turbulent diffusion term is 
anisotropic. Hence, the logic used in deriving (80) is virtually the same as that used in 

* MFI in the limit of two-dimensional turbulence can be satisfied identically by (76) if oio = — 3«9 + 12; 
see Speziale (1987a). 
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deriving the modeled ^-transport equation for the K — e model. Near solid boundaries, 
anisotropic corrections to (79) have been proposed that are typically of the algebraic form 
(see Hanjalic and Launder 1976) 

eij = ^eSij + 2ef t bij (81) 

where f a is a function of the turbulence Reynolds number Re f = K 2 / ve. Equation (81) - 
which can be thought of as a first-order Taylor expansion about a state of isotropic turbulence 
- is solved in conjunction with (80) where the model coefficients are taken to be functions of 
Re t as a solid boundary is approached (cf. Hanjalic and Launder 1976). As an alternative 
to (81), the isotropic form (79) can be used in a wall bounded flow if suitable wall functions 
are used to bridge the outer and inner flows. 

One major weakness of the models (80) - (81) is their neglect of rotational strains. For 
example, in a rotating isotropic turbulence, the modeled e-transport equation (80) yields the 
same decay rate independent of the rotation rate of the reference frame. In stark contrast 
to this result, physical and numerical experiments indicate that the decay rate of the tur- 
bulent kinetic energy can be considerably reduced by a system rotation - the inertial waves 
generated by the rotation disturb the phase coherence needed to cascade energy from the 
large scales to the small scales (see Wigeland and Nagib 1978 and Speziale, Mansour and 
Rogallo 1987). A variety of modifications to (80) have been proposed during the last decade 
to account for rotational strains (see Pope 1978, Hanjalic and Launder 1980, and Bardina, 
Ferziger and Rogallo 1985). However, these modifications have tended to be “one problem” 
corrections which gave rise to difficulties when other flows were considered. It was recently 
shown by the author that all of these modified e-transport equations are more ill-behaved 
than the standard model (80) for general homogeneous turbulent flows in a rotating frame 
(e.g., they fail to properly account for the stabilizing effect of a strong system rotation on a 
homogeneously strained turbulent flow; see Speziale 1989b). 

At this point it should be mentioned that in the second-order closure models of Mellor 
and co-workers, the dissipation rate is modeled as in equation (55) and a modeled transport 
equation for the integral length scale (70) is solved which is identical in form to (71). When 
this model has been applied to wall bounded turbulent flows it has typically been used in 
conjunction with wall functions. In addition, it should also be mentioned that second-order 
closure models along the lines of the K — w model of Wilcox and co-workers have been 
considered (here a modeled transport equation for the reciprocal time scale u> = e/K is 
solved; cf. Wilcox 1988). 

In order to complete these second-order closures, a model for the third-order diffusion 
correlation Cijk is needed. Since this is a third-order moment, the simplifying assumption 
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of gradient transport (which is generally valid only when there is a clear-cut separation of 
scales) is typically made. Hence, all of the commonly used second-order closures are based 
on models for C{jk of the form: 

n d T lm 

where the diffusion tensor Cijkimn can depend anisotropically on ry. For many incompressible 
turbulent flows, the pressure diffusion terms in Cijk can be neglected in comparison to the 
triple velocity correlation u t UjUk. Then, the symmetry of C^k under an interchange of any 
of its three indices immediately yields the form 





, drjk 
+ Tim 9x m 


“f* Tfcm 



(82) 


which was first obtained by Launder, Reece and Rodi (1975) via an alternative analysis based 
on the transport equation for v^UjUk ■ Equation (82) is sometimes used in its isotropized form 


r „ . *sA 

* jfe 3 * e \ dxi dxj dxk ) 


(83) 


(cf. Mellor and Herring 1973). The constant C, was chosen to be 0.11 by Launder, Reece 
and Rodi (1975) based on comparisons with experiments on thin shear flows. 


Examples 

Now, by the use of some illustrative examples, a case will be made for the superior 
predictive capabilities of second-order closures in comparison to zero, one and two equation 
models. First, to demonstrate the ability of second-order closure models to describe Reynolds 
stress relaxation effects, we will consider the return to isotropy problem. In this problem, 
an initially anisotropic, homogeneous turbulence - generated by the application of constant 
mean velocity gradients — gradually relaxes to a state of isotropy after the mean velocity 
gradients are removed. By introducing the transformed dimensionless time t (where dr = 
edt/2K), the modeled Reynolds stress transport equation can be written in the equivalent 
form 

^ + Aj (84) 

where Aij is the dimensionless slow pressure strain correlation. Since the rapid pressure- 
strain and transport terms vanish in this problem - and since the dissipation rate can be 
absorbed into the dimensionless time r - only a model for the slow pressure strain correlation 
is needed as indicated in (84). In Figure 5, the predictions of the LRR model (where 
Aij = —C\bij and the Rotta constant C\ = 3.0) for the time evolution of the second invariant 
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of the anistropy tensor are compared with the experimental data of Choi and Lumley (1984) 
for the relaxation from plane strain case. It is clear from this figure that this simple second- 
order closure model does a reasonably good job in reproducing the experimental trends 
which show a gradual return to isotropy where II — > 0 as t — * oo. This is in considerable 
contrast to eddy viscosity (or nonlinear algebraic stress) models which erroneously predict 
that II = 0 for r > 0! Further improvements can be obtained with second-order closures 
based on nonlinear models for the slow pressure strain correlation. A simple quadratic model 
was recently proposed by Sarkar and Speziale (1990) where 

Aj = -CAi + C,(Mw - i lib) (85) 

with Ci = 3.4 and C 2 = 4.2. This model does a better job in reproducing the trends of 
the Choi and Lumley (1984) experiment (see Figure 6). Most notably, the quadratic model 
(85) yields curved trajectories in the £ — r] phase space (where £ = = 7/5 ) that are 

well within the range of experimental data; any linear or quasilinear model where C 2 = 0 
erroneously yields straight line trajectories in the £-77 phase space as clearly shown in Figure 
6 . 

As alluded to earlier, second-order closure models perform far better than eddy viscosity 
models in rotating turbulent flows. To illustrate this point, a comparison of the predictions 
of the standard K — e model and the Launder, Reece and Rodi model will be made for the 
problem of homogeneous turbulent shear flow in a rotating frame. This problem represents 
a non-trivial test of turbulence models since a system rotation can have either a stabilizing 
or destabilizing effect on turbulent shear flow. The most basic type of plane shear flow in a 

s 0 \ 

0 0, fli = (0,0,0) (86) 

OO/ 

(see Figure 7). For the case of pure shear flow (0 = 0), the Launder, Reece and Rodi model 
yields substantially improved predictions over the K — e model for the equilibrium values 
of b{j and SKje as shown in Table 1. Since the standard K — £ model is frame- indifferent, 
it erroneously yields solutions for rotating shear flow that are independent of 0. Second- 
order closure models, on the other hand, yield rotationally-dependent solutions due to the 
effect of the Coriolis acceleration. For any homogeneous turbulent flow in a rotating frame, 
second-order closure models take the form (cf. Speziale 1989a) 


rotating frame will be discussed where 


duj 

dxj 


= 0 


$1% ’ &U • 

t \j = — T,fc J — Tjk~ b H*j — e *j — 2(Tj*e m *jf2 m + 7jfce m fcift m ) 

Ox k ox k 


(87) 
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where the mean vorticity tensor JJSij in the model for II*j (see equation (78)) is replaced with 
the intrinsic mean vorticity tensor Wij defined in (32). The equations of motion for the LRR 
model are obtained by substituting (86) into (87) and the modeled e-transport equation 


e = 


-C. i 


e dU{ 

K Tii 7hj 



( 88 ) 


which is not directly affected by rotations. A complete dynamical systems analysis of these 
nonlinear ordinary differential equations — which are typically solved for initial conditions 
that correspond to a state of isotropic turbulence - was conducted recently by Speziale 
and Mac Giolla Mhuiris 1989a. It was found that e/SK and have finite equilibrium 
values that are independent of the initial conditions and only depend on ft and S through 
the dimensionless ratio ft/S. There are two equilibrium solutions for (e/SK) <x> : one where 
( e/SK)oo = 0 which exists for all ft/5 and one where (e/SK) > 0 which only exists 
for an intermediate band of f l/S (see Figure 8(a)). The trivial equilibrium solution is 
predominantly associated with solutions for K and e that undergo a power law decay with 
time; the non-zero equilibrium solution (ellipse ACB on the bifurcation diagram shown in 
Figure 8(a)) is associated with unstable flow wherein K and e undergo an exponential time 
growth at the same rate. The two solutions exchange stabilities in the interval AB (i.e., 
this is a degenerate type of transcritical bifurcation). In stark contrast to these results, 
the standard K — e model erroneously predicts the same equilibrium value for (e/ S K)^ 
independent of Q/S (see Figure 8(b)). In Figures 9(a) - (c), the time evolution of the 
turbulent kinetic energy predicted by the standard K — e model and the LRR model are 
compared with the large-eddy simulations of Bardina, Ferziger and Reynolds (1983). It is 
clear that the second-order closure model is able to properly account for the stabilizing or 
destabilizing effect of rotations on shear flow whereas the K — e model erroneously predicts 
results that are independent of the rotation rate ft. The LRR model predicts that there 
is unstable flow (where K and e grow exponentially) only for rotation rates lying in the 
intermediate range — 0.1 £ ft/5 ^ 0.39 whereas linear stability analyses indicate unstable 
flow for 0 < ft /S < 0.5. Similar improved results using second-order closures have been 
recently obtained by Gatski and Savill (1989) for curved homogeneous shear flow. 

Finally, an example of an inhomogeneous wall-bounded turbulent flow will be given. The 
problem of rotating channel flow recently considered by Launder, Tselepidakis and Younis 
(1987) represents a challenging example. In this problem a turbulent channel flow is sub- 
jected to a steady spanwise rotation (see Figure 10). Physical and numerical experiments 
(see Johnston, Halleen and Lezius 1972 and Kim 1983) indicate that Coriolis forces arising 
from a system rotation cause the mean velocity profile u(y) to become asymmetric about the 
channel centerline. In Figure 11, the mean velocity profile computed by Launder, Tselepi- 
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dakis and Younis (1987) using the Gibson and Launder (1978) second-order closure model is 
compared with the results of the K —e model and the experimental data of Johnston, Halleen 
and Lezius (1972) for a Reynolds number Re = 11,500 and a rotation number Ro = 0.21. 
From this figure, it is clear that the second-order closure model yields a highly asymmetric 
mean velocity profile that is well within the range of the experimental data. The standard 
K — t model erroneously predicts the same symmetric mean velocity profile as in an inertial 
framing (where ft = 0) as shown in Figure 11. Comparable improvements in the prediction 
of curved turbulent shear flows have been obtained by Gibson and Rodi (1981) and Gibson 
and Younis (1986) using second-order closure models. Likewise, turbulent flows involving 
buoyancy effects have been shown to be better described by second-order closure models (cf. 
Zeman and Lumley 1976, 1979). In these problems, the Coriolis terms on the r.h.s. of the 
Reynolds stress transport equation (87) are replaced with the body force term 

Msl TVj + Si5x) (89) 

where j3 is the coefficient of thermal expansion and gi is the acceleration due to gravity. 
The temperature- velocity correlation T'u[ (also called the Reynolds heat flux) is modeled 
by a gradient transport hypothesis or is obtained from a modeled version of its transport 
equation. 

While second-order closure models constitute, by far, the most promising approach in 
Reynolds stress modeling, it must be said that they have not progressed to the point where 
reliable quantitative predictions can be made for a variety of turbulent flows. To illustrate 
this point we again cite the case of rotating shear flow. As shown earlier, the phase space 
portrait of second-order closures is far superior to that of any two-equation model of the 
eddy viscosity type (i.e., the second-order closures properly predict that there is unstable 
flow only for an intermediate band of rotation rates; see Figure 8). However, the specific 
quantitative predictions of a wide variety of existing second-order closures were recently 
shown by Speziale, Gatski and Mac Giolla Mhuiris (1989) to be highly contradictory in 
rotating shear flow for a significant range of f l/S (see Figure 12). Comparable problems 
with the reliability of predictions when second-order closure models are integrated directly 
to a solid boundary persist so that a variety of modifications - which usually involve the 
introduction of empirical wall damping that depends on the turbulence Reynolds number 
as well as the unit normal to the wall - continue to be proposed along alternative lines (cf. 
Launder and Shima 1989, Lai and So 1989 and Shih and Mansour 1990). 

In the opinion of the author, there are two major areas of development that are direly 
needed in order to improve the predictive capabilities of second-order closures: 

(1) The introduction of improved transport models for the turbulence length scale which 
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incorporate at least some limited two-point and directional information (e.g., through 
some appropriate integral of the two-point velocity correlation tensor Rij). In conjunc- 
tion with this research, the use of gradient transport models should be re-examined. 
Although Donaldson and Sandri (1981) developed a tensor length scale along these 
lines, it was recently shown by Speziale (1989b) that the specific form of the model 
that they chose can be collapsed to the standard e-transport model in homogeneous 
flows. 

(2) The need for asymptotically consistent low turbulence Reynolds number extensions of 
existing models that can be robustly integrated to a solid boundary. Existing models 
use ad hoc damping functions based on Re t and have an implicit dependence on the 
unit normal to the wall which does not allow for the proper treatment of geometrical 
discontinuities such as those that occur in the square duct or back-step problems. 
Furthermore, the nonlinear effect of both rotational and irrotational strains need to 
be accounted for in the modeling of near wall anisotropies in the dissipation. 


CONCLUDING REMARKS 

There has been a tendency to be overly pessimistic about the progress made in Reynolds 
stress modeling during the past few decades. It must be remembered that the first complete 
Reynolds stress models - cast in tensor form and supplemented only with initial and bound- 
ary conditions - were developed less than twenty years ago. Progress was at first stymied by 
the lack of adequate computational power to properly explore full Reynolds stress closures in 
non-trivial turbulent flows - a deficiency that was not overcome until the late 1960’s. Then, 
by 1980 - with an enormous increase in computer capacity - efforts were shifted toward 
direct and large eddy simulations of the Navier-Stokes equations. Furthermore, the interest 
in coherent structures (cf. Hussain 1983) and alternative theoretical approaches based on 
nonlinear dynamics (e.g., period doubling bifurcations as a route to chaos; cf. Swinney and 
Gollub 1981) that crystalized during the late 1970’s have also shifted attention away from 
Reynolds stress modeling, as well as the general statistical approach for that matter. While 
progress has been slow, this is due in large measure to how intrinsically difficult the problem 
is. The fact that real progress has been made, however, cannot be denied. Many of the 
turbulent flows considered in the last section - which were solved without the introduction 
of any further ad hoc empiricisms - could not be properly analyzed by the Reynolds stress 
models that were available before 1970. 

Some discussion is warranted concerning the goals and limitations of Reynolds stress mod- 
eling. Under the best of circumstances, Reynolds stress models can only provide accurate in- 
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formation about first and second one-point moments (e.g., the mean velocity, mean pressure, 
and turbulence intensity) which, quite often, is all that is needed for design purposes. Since 
Reynolds stress modeling constitutes a low-order one-point closure, it intrinsically cannot 
provide detailed information about flow structures. Furthermore, since spectral information 
needs to be indirectly built into Reynolds stress models, a given model cannot be expected 
to perform well in a variety of turbulent flows where the spectrum of the energy containing 
eddies is changing dramatically. However, to criticize Reynolds stress models purely on the 
grounds that they are not based rigorously on solutions of the full Navier- Stokes equations 
would be as childish as criticizing exact solutions of the Navier-Stokes equations for not 
being rigorously derivable from the Boltzmann equation or, for that matter, from quantum 
mechanics. The more appropriate question is whether or not a Reynolds stress model can 
be developed that will provide adequate engineering answers for the mean velocity, mean 
pressure and turbulence intensities in a significant range of turbulent flows that are of techno- 
logical interest. To obtain accurate predictions for these quantities in all possible turbulent 
flows will probably require nothing short of solving the full Navier-Stokes equations. Such a 
task will not be achievable in the foreseeable future, if ever at all (cf. Hussaini, Speziale and 
Zang 1989). To gain an appreciation for the magnitude of such an endeavor, consider the 
fact that economically feasible direct simulations of turbulent pipe flow at a Reynolds num- 
ber of 500,000 - a turbulent flow which, although non-trivial, is far from the most difficult 
encountered - would require a computer that is 10 million times faster than the Cray YMP! 

While second-order closures represent the most promising approach in Reynolds stress 
modeling, much work remains to be done. The two problem areas mentioned in the previous 
section - namely, the development of transport models for an anisotropic integral length scale 
and the development of more asymptotically consistent methods for the integration of second- 
order closures to a solid boundary - are of utmost importance. In fact, the latter issue of near 
wall modeling is so crucial that deficiencies in it - along with associated numerical stiffness 
problems - are primarily responsible for the somewhat misleading critical evaluations of 
second-order closures that arose out of the 1980-81 AFOSR-HTTM Stanford Conference on 
Turbulence (see Kline, Cantwell and Lilley 1981). Another area that urgently needs attention 
is the second-order closure modeling of compressible turbulent flows. Until recently, most 
compressible second-order closure modeling has consisted of Favre-averaged, variable-density 
extensions of the incompressible models (cf. Cebeci and Smith 1974). However, with the 
current thrust in compressible second-order modeling at NASA Langley and NASA Ames, 
some new compressible modeling ideas - such as dilatational dissipation - have come to the 
forefront (see Sarkar et al. 1989 and Zeman 1990). Much more work in this area is needed, 
however. 
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Reynolds stress modeling should continue to steadily progress - complementing numerical 
simulations of the Navier-Stokes equations and alternative theoretical approaches. In fact, 
with anticipated improvements in computer capacity, direct numerical simulations should 
begin to play a pivotal role in the screening and calibration of turbulence models. Further- 
more, from the theoretical side, statistical mechanics approaches such as RNG could be of 
considerable future use in the formulation of new models (unfortunately, at their current 
stage of development, it does not appear that they can reliably calibrate turbulence models 
for use in complex flows). Although Reynolds stress models provide information only about 
a limited facet of turbulence, this information can have such important scientific and engi- 
neering applications that they are likely to remain a part of turbulence research for many 
years to come. 
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Equilibrium 

Values 

S t an dard 
K—e Model 

Launder, Reece 
& Rodi Model 

Experiments 

(6 n)oo 

0 

0.193 

0.201 

( 622)00 

0 

-0.096 

-0.147 

( 612)00 

-0.217 

-0.185 

-0.150 

(SK/e) oo 

4.82 

5.65 

6.08 


Table 1. Comparison of the predictions of the standard K — e model and the Launder, 
Reece, and Rodi model with the experiments of Tavoularis and Corrsin (1981) 
on homogeneous shear flow. 


43 








46 


Figure 2. Comparison of the predictions of the standard K — € model with the ex- 

periments of Kim, Kline, and Johnston (1980) for turbulent flow past a 
backward facing step: (a) mean flow streamlines, (b) turbulence intensity 
profiles. 
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turbulence intensity profiles. 
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Figure 4. Comparison of the predictions of the nonlinear K — £ model of Speziale 

(1987b) with the experiments of Kim, Kline, and Johnston (1980) for tur- 
bulent flow past a backward facing step: (a) mean flow streamlines, (b) 
turbulence intensity profiles. 
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Figure 5. Time evolution of the second invariant of the anisotropy tensor for the 

return to isotropy problem: Comparison of the predictions of the Launder, 
Reece, and Rodi model with the experimental data of Choi and Lumley 
(1984) for the relaxation from plane strain. 


50 



£=///'/3 









(a) K-t Model 


b LRR Model 


n I/S = 0, 0.25 4 -0.50 


Time evolution of the turbulent kinetic energy in rotating shear flow for 
£o/SK 0 = 0.296: (a) standard K — e model, (b) Launder, Reece, and Rodi 
model, and (c) large-eddy simulations of Bardina, Ferziger, and Reynolds 







Figure 10. Fully- developed turbulent channel flow in a rotating frame. 
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Figure 11. 



y/D 


Comparison of the mean velocity predictions of the second-order closure 
model of Gibson and Launder (1978) and the standard K — t model with 
the experimental data of Johnston, Halleen, and Lezius (1972) on rotating 
channel flow. 
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Figure 12. Comparison of the predictions of a variety of second-order closure models 

for the time evolution of the turbulent kinetic energy in rotating shear flow; 
0,/S = 0.25, Cq/SKq = 0.296. LES = large-eddy simulations of Bardina, 
Ferziger, and Reynolds (1983); LRR = Launder, Reece, and Rodi model; 
RK = Rotta Kolmogorov model of Mellor and Herring (1973); FLT = Fu, 
Launder, and Tselepidakis (1987) model; RNG = renormalization group 
model of Yakhot and Orszag; SL = Shih and Lumley (1985) model. 
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